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Abstract 

Bulk nuclear observables such as charge radii and binding energies are well described by both 
nonrelativistic and covariant mean-field models. However, predictions of neutron radii, which are 
not tightly constrained by reliable data, vary signihcantly. The nature of this variation is investi¬ 
gated using correlations between basic properties of the models and the neutron skin thickness in 
lead. The results suggest that conventional covariant models are too limited. The study is guided 
by principles and insights of effective field theory (EFT), such as power counting, and the relation 
of mean-field models to a more general EFT approach to nuclei is discussed. 
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I. INTRODUCTION 


Bulk nuclear observables such as charge radii and binding energies are well described by 
both nonrelativistic and covariant mean-held models ||I|, Yet the “best ht” calculations 
of each type in the literature disagree substantially on neutron radii The spread in 

mean-held predictions and a recent analysis of experimental determinations suggest that 
neutron radii in medium to heavy nuclei are inadequately determined at present by either 
theory or experiment. Better values for neutron radii could have widespread impact, from 
reducing uncertainties in atomic parity violation experiments p] to constraining properties 
of the surface crust of neutron stars |^. More generally, the lack of control over neutron 
radius predictions implies large uncertainties in extrapolations to neutron- or proton-rich 
nuclei, which play an important role in nuclear astrophysics |^. 

There are arguments in the literature that mean-held models of nuclei are inadequate 
because the mean-held approximation is not valid at ordinary nuclear densities or that essen¬ 
tial correlations are excluded [^]. However, these are misleading arguments, which implicitly 
assume that one is performing a Hartree or Hartree-Fock calculation of some underlying 
interaction. In fact, it is more appropriate to identify these calculations with Kohn-Sham 
density functional theory (DFT) 0, ^ P, in which the full Hartree plus exchange-correlation 
functionals (and not Hartree alone!) are approximated by a parametrized form.* The im¬ 
portant observation is that mean-held models are specihed by a universal energy functional 
of the density, which is minimized iteratively for each nucleus through the introduction of 
auxiliary single-particle orbitals. This is the structure of a Kohn-Sham DFT. The solution 
procedure takes a small fraction of the computational ehort of conventional direct many- 
body solutions. The universality is important because the difficulties in hnding solutions 
to many diherent nuclei are shifted to constructing the functional once; if the functional is 
easily evaluated then the cost of subsequent applications to each additional hnite nucleus is 
small, and one can scale to large numbers of nucleons. 

“Mean held” in this context really implies a limited form of the analytic and nonlocal 
structure that can appear in the exact energy functional. We interpret the mean-held mod¬ 
els as a form of Kohn-Sham DFT in which the functional is approximated by a type of 
generalized gradient expansion |0,113. From this point of view, the underlying framework 


is completely general and the limitations of present-day mean-held calculations could be 
systematically removed without losing the computational advantages of the approach. Defi¬ 
ciencies in the mean-field model predictions should therefore be traced to deficiencies in the 
approximate functional rather than to a breakdown of the underlying framework. 

In this paper, we adapt this philosophy to a study of mean-held neutron radii. The nature 
of the variation in radii is investigated using correlations between basic characteristics of the 
models and the predicted neutron skin thickness in lead. In this way, we can identify 
how properties of the mean-held energy functionals are related to observables and identify 
possible dehciencies in the functional. The study is guided by the principles and insights 
of ehective held theory (EFT), with an eye on the ultimate goal of embedding mean-held 
models into a more general DFT/EFT approach to nuclei |^. With an EFT approach 
one can exploit power counting to systematically approximate the energy functional. In the 
present context, power counting means that one can estimate the contribution to the energy 


The observation that mean-held models can be interpreted as Kohn-Sham functionals has been made by 
Brack for Skyrme models 0 and by Schmid et al. for covariant models j0 0. 
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of individual terms in the energy functional, even those excluded by the model truncation. 
The promise of this approach is supported by empirical evidence that EFT power counting 
in mean-held models is robust 130,00. 

In principle, a detailed analysis during the fitting process would best reveal how observ¬ 
ables such as the neutron radius are correlated with other properties of the model. As an 
alternative, we consider a large set of models all ht to roughly the same data at about the 
same accuracy, and then examine the correlations graphically. This strategy has been used 
with covariant models to study deformations in nuclei [M and spin-orbit splittings pi 


Brown has used such correlations with the Skyrme model recently to address many of the 
same issues about neutron radii and Typel and Brown have extended these results to 


conventional covariant models In our analysis, we use conventional parameter sets from 
the literature but also generalized functionals developed for EFT studies. The latter provide 
a wide class of models that provide very good hts to the usual bulk nuclear data but which 
also span a wide range of observables that are not well constrained by this data. 

We hnd that the neutron skin thickness is dominantly correlated with the bulk symmetry 
energy and its density dependence. Thus the properties of asymmetric inhnite nuclear 
matter (as opposed to surface properties) are most important. Our conclusions on the 
physics determining the skin are in accord with many of the observations in Ref. , which 


are based on applying a macroscopic model and comparing to two particular mean-field 
models, and by Brown and Brown and Typel (Note: many of the same issues were 
discussed long ago by Bodmer who used a semi-empirical statistical method |^. Other 
recent studies of the neutron skin in mean-held models include Refs. 


M 1^, 1^, p9 


The differences between predictions from standard Skyrme and covariant models can then 
be traced to differences in the density dependence of the symmetry energy. A comparison of 
the symmetry energy derived from more microscopic approaches that are tied to free-space 
nucleon-nucleon scattering implies that conventional covariant models based on heavy meson 
interactions are too simple. 

The plan of the paper is as follows. In Sect. ^ we describe the mean-held functionals and 
and in Sect. JTT|the EFT power counting. The correlation analysis is carried out in Sect. |V 
and the consequences discussed in Sect. ^ Sect. ^ 


IS a summary. 


II. OVERVIEW OF MEAN-FIELD ENERGY FUNCTIONALS 

There are two major classes of mean-held models in general use, which can be formulated 
in terms of nonrelativistic (Skyrme) or covariant energy functionals. The covariant models 
can be further subdivided into meson and point-coupling models. We will compare all 
three types in our correlation analysis. These models and their implementations are well 

|36| , ^ ^ |40| and we 

will only discuss here details relevant to the present investigation. We focus directly on 
energy functionals instead of Hamiltonians or Lagrangians to stress the connection to a 
more general formulation in terms of density functional theory. The usual construction of 
the functional starts by postulating interactions (either a Skyrme ehective interaction or a 
Lagrangian) and then treats them in the Hartree approximation.! 


documented in the literature i, i 0, 0, 


34, 35 


! Note that the Skyrme models actually apply a Hartree-Fock approximation, but since the forces are zero 
range, this doesn’t affect the form of terms in the energy functional. 
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The Skyrme models are based on nonrelativistic zero-range, density-dependent interac¬ 
tions. The energy functional most commonly used in the present-day literature has evolved 
only slightly from the original form proposed; an explicit expression for the interaction is 
given in Ref. [Q]. The corresponding functional can be written as powers and gradients of 
the isoscalar density p = Pp + Pn and isovector density p^ = Pp — pn, the corresponding 
kinetic energy densities r and and spin-orbit currents. As such it can be considered a 
Kohn-Sham functional constructed as a generalized gradient expansion |^, |^. Of partic¬ 
ular note are terms proportional to and p“p 3 , which in modern Skyrme parameter sets 
appear with fractional a. Although extended Skyrme models have been considered PT| , 
the density expansion stops with these terms. There are usually four isovector parameters 
(which are linear combinations of the traditional parameters xq-x^, and Wq [| 1 |), cor¬ 


responding to P 3 T 3 , and (Vp 3 )^ terms (there are also isovector terms involving the 

spin-orbit currents). While the Skyrme interaction has been motivated in terms of underly¬ 
ing nucleon-nucleon interactions |^, the connection has only been made qualitatively 
and the parameters are treated phenomenologically. A goal of the EFT/DFT program 
advocated here will be to revisit this connection and seek systematic improvement. 

Many Skyrme parameter sets are available in the literature 


22, 32 


I, representing 

a wide range of fits to bulk nuclear properties, although the fits are not all to the same set of 
observables and are not of equal quality. The Skyrme parameter set developed by Brown |I| 
is fit to the widest range of properties and serves as the quality-of-fit benchmark for this class 
of models. We will use most of the sets from the literature plus many new ones generated for 
this investigation. The new sets are fit to binding energies, charge densities, and spin-orbit 
splittings of doubly magic nuclei as described for covariant models in Refs. |T^ and [I^ . 


However, a value for the neutron radius in lead was also specified so as to generate sets 
that cover as large a range in neutron radii as is compatible with a good £t to the other 
properties. Only new sets with a quality of fit comparable to Ref. were used in our 
analysis. 

The most commonly used covariant mean-held models are derived from Lagrangians 
with nucleons coupled to heavy mesons and are characterized by large, isoscalar, mean 
meson helds that are several hundred MeV in magnitude 


35 


Meson exchange 

in these models is spacelike and far off-shell, so the corresponding parameters are at best 
only roughly related to the physical meson spectrum. We describe these energy functionals 
as “covariant” rather than “relativistic,” as is conventional in the literature, since the latter 
suggests the importance of kinematic relativity—large velocities and large p/M —whereas 
these corrections are small for nuclei. Rather it is the potentials where relativity is most 
important, and this is manifested in a Lorentz covariant representation, which maintains the 
distinction between Lorentz scalars and zeroth components of four-vectors (It is not 

Lorentz invariant, since the four velocity of a nucleus or nuclear matter dehnes a preferred 
frame.) 

In practice, the form of the energy functional has been limited to terms that can arise 
from the one-loop or Hartree approximation to a Lorentz invariant Lagrangian. In this sense 
the functionals are less general than the Skyrme functionals. This comparison is somewhat 
misleading, however, since a functional built from integral powers of both Lorentz scalar 
and vector densities will, when reduced to nonrelativistic form, contain terms non-analytic 
in the density (i.e., analogous to the term for particular values of a) [^, 0|. There are 
several standard covariant-model parameter sets in the literature |^, and many more 
have been generated recently in the course of EFT investigations [^, We use them 
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all in the present investigation. In conventional mean-field parametrizations, the model 
is usually truncated to minimize the parameters needed to provide a good £t, with the 
hope of maximizing predictive power. Terms through the fourth power of the scalar mean 
field are included but only through the second power of the other mean fields; there is 
only one non-gradient isovector parameter. In an EFT-inspired analysis, one uses power 
counting arguments to organize the energy functional and reveal how many parameters can 
be fixed from the data. The EFT parameter sets include more interaction terms but are 
still not complete. For example, a meson corresponding to isovector scalar exchange was not 
considered because its effects were found to be small in previous investigations. In Ref. [TO , 
this is argued to have an important effect on the symmetry energy. 

An alternative to explicit mesons are “point coupling” (PC) models, which feature contact 
interactions and derivatives among nucleons only rather than meson exchange [O, M ^ . 


The energy functional here is composed of powers and gradients of Lorenz scalar (p^), vector 

H 


(Pb) and tensor densities and currents 


Point coupling models are somewhat more flex¬ 


ible than covariant meson models because the assumption of an underlying meson exchange 
does not constrain terms in the energy functional. Thus a term such as (Vps)^ would have 
a dehnite sign if based on an underlying meson exchange but could have either sign in a 
point-coupling model. 

In both types of covariant models, the isovector interactions in conventional formulations 
are quite limited, with the form of simple rho meson exchange or its point-coupling equiv¬ 
alent. Since the corresponding mass is largely unconstrained by hts to nuclear properties, 
there is really only one free parameter. This is in contrast to conventional Skyrme interac¬ 
tions, which have up to four parameters controlling isovector contributions (although two 
are largely irrelevant for the properties discussed here |M). 


III. POWER COUNTING IN THE ENERGY FUNCTIONAL 


An effective field theory (EFT) describes low-energy physics with low-energy degrees of 
freedom. Underlying EFT is the principle that while the short-distance, ultraviolet behavior 
of the effective theory may be incorrect, it can be corrected systematically by the renormal¬ 
ization of local operators (“counterterms”) There is considerable freedom in the choice 
of representation of these operators; while one choice is not more “correct” than another, 
it may be more convergent or easier to calculate. It is in this context that we compare 
nonrelativistic and covariant energy functionals of the density for nuclei. A summary of the 
EFT philosophy as applied to mean-held models is given in Ref. [^. For our purposes, the 
most important concept is that of power counting in the energy functional, which associates 
a natural size to every term based on conjectured underlying scales. 

Will an EFT approach be useful for heavy nuclei? An EFT relies on the separation of 
scales, so we should identify the momentum scale A that divides short from long-distance 
physics. Long-distance physics must be treated explicitly while short-distance physics can 
be reproduced with local operators. For the latter we might expect a density expansion in 
powers of k-p/A, where kp is the Fermi momentum. (This is manifest in the treatment of 
a dilute Fermi system but has not been shown for nuclear systems.) There are two 
immediate problems. First, nuclear matter equilibrium kp is signihcantly larger than the 
pion mass, which implies that an explicit treatment of pion physics as long-range physics is 
essential for a useful EFT. Second, nuclear matter saturation might only occur when kp ~ A, 
because only then could different terms in the density expansion compensate each other. 
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However, the bulk physics of heavy nuclei is dominated by isoscalar physics (see below), 
and the dominant pion contributions appear to be two-pion physics, which is associated 

The success in using a exchange in two-nucleon 


36 


with the Lorentz scalar “a” meson 
potentials implies that the low-energy tail of the two-pion contribution is not well resolved, 
particularly in nuclei.^ The consequence is that in an EFT, nuclear structure is dominated 
by short-distance contact terms (“low-energy constants”), which implies that a model with 
A 600 MeV without explicit pions is actually a good starting point. This argument also 
implies that isovector observables such as the neutron skin, which could be sensitive to 
one-pion-exchange physics, may be inadequately described. 

The potential problem with applying EFT to saturating systems, in which saturation 
is only possible at the breakdown scale, is apparently not realized in nuclei. In fact, by 
conventional analysis, equilibrium nuclear matter is an anomalously dilute Fermi system. 
One such analysis, by Jackson is based on saturation driven by a strongly repulsive 

short-range interaction (which we’ll call a “hard core”) with range c (so that A 1/c). The 
size i? of a saturating system scales with the number of particles A to the one-third power, 
which dehnes a characteristic length scale tq: 


R ~ roA^/3 . 


( 1 ) 


Jackson estimates the limits of saturation to be bounded by 0.552c < tq < 2.4c ^I|. This 
estimate is consistent with liquid ^He, where tq ~ 0.96c, but with a nuclear hard core at 
c = 0.4 fm, rg is actually beyond the bound: tq ^ 2.75c, which implies that saturation in 
nuclei is unconventional. One possible explanation for saturation in an EFT with a density 
expansion is that two low orders might largely cancel because of an anomalously small 
coefficient multiplying a leading term, while all higher orders follow the expansion in /cp/A. 
This is the mechanism for saturation in covariant mean-held models 


15 


We assume that a counting consistent with low-energy QCD, including pions as light 
degrees of freedom, is appropriate for our mean-held energy functionals. Georgi-Manohar 
naive dimensional analysis, or NDA, assigns appropriate powers of the pion decay constant 
Ri 93 MeV and a scale of non-goldstone boson physics A (about 600 MeV) to terms in 
a low-energy ehective Lagrangian of QCD |^, |5^. The counting for a generic term in a 
covariant meson Lagrangian is found from 


/3 lAW 



1 

1 

fgv^ 

V A3 J ' 

V A / m! 

VAJ nV 

1 A ) 


( 2 ) 


where ~ A//^ and Dirac and isospin operators are not shown. (See Ref. [T^ for a discussion 
of power counting with covariant meson Lagrangians including pions.) The assumption here 
is that this power counting, intended for terms in a Lagrangian, can be applied directly to the 
energy functional through the associations ps <->■ ip'ilj, Pe ^ V’V ^md so on. This connection 
is immediate if the mean-held model is viewed as being only a one-loop calculation but 
has not been fully clarihed in the more general context of DFT. Skyrme functionals follow 
an analogous power counting behavior with powers of p and r counted using p 
r and so on (see Ref. [|T^ for details). 


If we have identihed the appropriate underlying length scales, naturalness implies that the 
dimensionless constant (3 is of order unity. Given this result, we can estimate contributions 


t Indeed, one must work very carefully to discern the long-range chiral two-pion-exchange tail in the detailed 
two-nucleon scattering data . 
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FIG. 1: Contributions to the energy per particle in for two covariant point-coupling models 

from Ref. 0- Absolute values are shown and the filled symbols are net values for the sum of 
scalar and vector terms. The small symbols indicate estimates based on NDA [Eq. @], with the 
error bars corresponding to natural coefficients from 1/2 to 2. 



to the energy from any term in the mean-field energy fnnctional. For example, a local- 
density estimate (per nucleon) for a general isoscalar term with dimensionless coefficient (3 
can be made according to 




d?x(3 {ps)\p^y ^ /3A 


pa ) 


(3) 


where p/ and pe are naturally scaled scalar and baryon densities (e.g., p/ = ps/pA) (ps) 
is the average density in the nucleus being considered. We take l/2</9<2asa reasonable 
range of natural coefficients. Other local density prescriptions are detailed in Refs. [^] and 
pl| . Figure]^ shows these natural estimates with the associated error bars from the variation 
in (3 as small squares and triangles for Actual contributions to the energy per particle 

in lead for two covariant point-coupling models that £t nuclear properties are obtained by 
applying the Hellmann-Feynman theorem [M. The NDA estimates are validated by these 
results and many other studies 0 H, n, p, n, n. 

We will take as a guiding principle that NDA and naturalness are valid for mean-held 
energy functionals. Indeed, all of the best-fit models exhibit naturalness in all or almost all 
coefficients. A comparison of goodness-of-ht with this counting at different levels of trunca¬ 
tion reveals the limit of the energy in constraining terms. In particular, contributions of order 
1 MeV are barely constrained, so contributions from beyond the fourth power of isoscalar 
densities or meson mean helds are not distinguished from numerical noise [|1^ . These results 
explain why the conventional covariant mean-held truncation is phenomenologically success¬ 
ful. In Ref. [0, a complete analysis like this concluded that only six constants in covariant 
energy functionals are well determined by hts to the usual bulk properties of nuclei. These 
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includes four isoscalar terms, one gradient term, and one isovector term. Brown concludes 
that a similar number of parameters are determined in fits to Skyrme models that use only 
the standard observables []^. 

The conventional set of observables used to £t parameters directly constrains the isovector 
contribution only through the binding energy (i.e., the symmetry energy). Isovector con¬ 
tributions can be reliably estimated using {p^} ^ {Z — N)/A x (pb)- The leading isovector 
contribution in Fig. is therefore suppressed relative to the leading isoscalar contribution 
by the factor [{Z — N)/2A]‘^ (the extra 1/2 comes from |r 3 [^), which is small even in lead. 
This puts the leading term on the order of 1-2 MeV, which is just resolved by the hts (pg 
isoscalar contributions are not resolved). When more than one parameter is available, as 
shown in the hgure, only a linear combination is determined, and the subleading coefficient is 
borderline unnatural |^. While individual values for these couplings vary widely between 
different parameter sets, the linear combination is well determined. The power counting 
analysis supports previous observations in meson models that, at the mean-field level, the 


role of the p meson is simply to adjust the symmetry energy [^, 


We emphasize that despite the phenomenological success of the power counting, we do 
not as yet have a systematic EFT expansion for the energy functionals. In particular, 
long-range loop contributions may lead to essential nonlocalities and nonanalyticities in 
the energy functional. An underlying assumption that we are testing is that long-distance 
contributions do not disturb the power counting. One should not assume, however, that 
long distance degrees of freedom imply that the functional cannot take the form of powers of 
densities and gradients. For example, we note that Coulomb systems, which depend entirely 
on long range physics, are well described by energy functionals based on local density and 
gradient expansions H. 


IV. CORRELATION ANALYSIS 


Most observables in finite nuclei have a highly nonlinear and correlated dependence on 
the parameters of the energy functional. Predicting the correlations between features of 
the mean-held functional and observables such as the neutron radius or the neutron skin 
thickness is prone to error and plausible explanations often turn out to be incorrect (as 
illustrated below for hnite range effects). This is where power counting and EFT arguments 
can be illuminating, even in the absence of a systematic EFT expansion. One can also go 
wrong by considering only a few “best £t” parameter sets, which can lead to false conclusions 
because of the limited sampling. 

Our strategy is to use many parameter sets, each a good £t to bulk properties of doubly 
magic nuclei. We consider generalized versions of each type of mean-held model and look at 
correlations between observables and characteristics of the models. In order to hll in gaps 
that can obscure trends, we also force particular neutron radii, being careful to keep only 
sets with acceptably good hts to the standard properties. We leave open the question of 
whether other hnite nucleus observables (such as inelastic scattering to collective states) can 
further restrict the acceptable sets. 

In the correlation plots that follow, the shape of a symbol indicates the type of mean- 
held model: circles for Skyrme, squares for covariant meson, and triangles for covariant 
point coupling. The shading indicates the origin of the model. Black hlled symbols are 
standard parameter sets from the literature |^, I 


33, 34, 38 


], Grey hlled symbols 
are parameter sets generated for EFT investigations [O, M and diagonal striped symbols 






























are unpublished generalized EFT sets from Ref. [45|. Black and white checkered symbols are 


new parameter sets generated for this investigation. Other than this (partial) identihcation 
we treat all sets democratically. Note that not all models are plotted on all figures. 

Comparisons between neutron radii or neutron skins from any given Skyrme model and 
any given covariant model over a wide range of nucleon number A show that the differences 
are very systematic. Fluctuations in the difference for particular nuclei, which are affected 
by details such as the treatment of pairing in open-shell nuclei, are small compared to the 
difference itself We conclude that understanding the skin thickness for one convenient 

nucleus is sufficient at the resolution of this investigation. We focus on ^ospp gj^ce it has 
the largest skin of the doubly magic nuclei. 

The skin thickness is a more robust isovector observable than the neutron radius alone. 
In principle, the charge radius of Pb^°® is very well determined as are the factors that 
“remove” the charge form factor of the proton (and other smaller corrections). Thus, the 
point proton radius should be very well hxed, and therefore the skin thickness should 


be closely correlated to r„. In practice, when 


- Ur 


is plotted as a function of for 


conventional models, there is a considerable spread, as seen in Fig. g. This spread reflects 
slightly different form factor corrections used in the models and differences in how well 
(r^)ch in jg reproduced by the hts. We want to isolate isovector properties, so the 

skin thickness is more informative than the individual radii. Since the neutron radii are not 
included in the hts determining these sets, the scatter is primarily driven by the binding 
energies (which reflect the bulk symmetry energy). 

Thus, we focus on the neutron skin thickness, (r„ — r^), in ^ospp g^g g representative 
isovector “observable” for our study. We put observable in quotes here because, in fact, 
neither the point proton radius, the point neutron radius, nor their difference is actually an 
experimental observable. Although it is standard practice, the “unfolding” of a point proton 
radius from a charge radius measurement is necessarily model dependent. [See Ref. for 
a discussion from the EFT point of view of the ambiguities in other quantities (occupation 
numbers) that are often treated as observables.] We will assume that the natural size 
of this ambiguity is small compared to the size of the difference in neutron radii, but this 
should be investigated further. A consistent EFT framework will obviate the need to discuss 
point data, as illustrated in Ref. [T^, where nucleon charge form factors are built into the 


Lagrangian and energy functional. 

We also note that the neutron radius is very highly correlated with the neutron form factor 
at low momentum transfer. This correlation is important for a proposed parity violation 
experiment at Jefferson Lab, which seeks to make a hve percent measurement of the form 
factor at a momentum transfer of 0.45 fm“^. An analysis of the experiment is given by 
Horowitz et al. P|. The relationship between the form factor at this momentum transfer and 
the neutron radius in ^osp]^ jg gPown in Fig. |^; it is evident that the proposed measurement 
will directly constrain the neutron radius. 

These hgures already show that conventional covariant meson models (black squares) in 
general predict larger skin thicknesses than conventional Skyrme models (black circles). A 
plausible explanation for the larger neutron skin thickness in meson models compared to 
Skyrme is the hnite range of the interaction in the former. The argument is that hnite range 
meson exchange enables the neutron density to extend further than the proton density and 
still feel the attraction from the protons. However, this picture is not consistent with an EFT 
viewpoint. The resolution scale associated with meson exchange in the covariant models is 
A 600 MeV. For momenta small compared to this scale, the details of the meson exchange 
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FIG. 2; Neutron skin thickness vs. neutron radius for a wide variety of mean-field models, as 
described in the text. The shading indicates the origin of the model. Black filled symbols are 
standard parameter sets from the literature |^, 


, 1^, ^|. Grey filled symbols 

are parameter sets generated for EFT investigations [IT, and diagonal striped symbols are 
unpublished generalized EFT sets from Ref. [45|. 
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FIG. 3: The calculated neutron form factor (as defined in Ref. [^) for ^*^®Pb at momentum transfer 
Q = 0.45 fm“^ vs. the point neutron radius in ^^®Pb for a wide variety of Skyrme and covariant 
mean-field models, as described in the text and in the caption to Fig. ||. 
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and the substructure of the hadrons are not resolved. As a result, this short-range physics 
can be incorporated into the coefficients of operators organized as a derivative expansion. 
The consequence for low-momentum properties such as the rms radii is that one cannot 
distinguish interactions with explicit ranges of order from interactions composed of S- 
functions plus derivatives. In principle, one needs an infinite number of terms for complete 
equivalence; in practice power counting tells us we can truncate sharply. 

At tree level, which corresponds to treating the functional as arising from the Hartree 
approximation, the analysis is trivial. The static propagator of a meson with mass m ~ A 
can be expanded in a Taylor series about |qp = = 0: 


m 


q2 _|_ ^2 


= 1 - 




+ 0{q^/m^) 1 — 




( 4 ) 


In general, there are two ingredients to the conclusion that the higher gradient terms can 
be neglected in an energy functional: i) -C A^ for typical and ii) the coefficient is of 
order unity [as in Eq. (^]. A predictive power counting prescription relies on both. For 
nuclei such as lead, (g^) is small in the interior and is largest in the surface, where the scale 
is set by the surface thickness a. Typically l/aA 1/5 JT^. The power counting analysis 
in Ref. [16| shows that resolving individual mass terms, which would mean determining two 


terms in each isoscalar meson propagator, is not possible (that is, only one gradient term is 
constrained by the £t). 

We can test this conclusion directly by varying the heavy meson masses. In Fig. the 
skin thicknesses of three nuclei are plotted against the isoscalar vector meson mass in a 
covariant meson model. As the vector mass was varied over a factor of two, the model 
parameters were adjusted to achieve roughly the same goodness of £t. The corresponding 
scalar mass is shown on the top axis (note that it varies much less than a factor of two). 
While strict attention to the of the £t favors a vector meson mass between 700 and 
800 MeV, the variation over the entire range is remarkably small. This conhrms that the 
range of the individual mesons is not well determined by the fit to nuclei and that a gradient 
expansion is more than adequate for this observable. 

The concept of optimal parameters for covariant models was introduced in Refs. 0,0 

m 


and further exploited in Ref. 


This is a reorganization of the power counting hierarchy 


by introducing linear combinations of parameters based on the close cancellations in covari¬ 
ant models between scalar and vector contributions. By rewriting in terms of the optimal 
parameters, the number of parameters actually determined by the data is manifested. For 
example, the combination: 


m? m?. 


( 5 ) 


is much better determined than each of these terms individually. The choice of optimal 
parameters is not unique and will depend on the quantity to be studied. To identify appro¬ 
priate optimal parameters for studying masses (gradients) in covariant meson models, we 
note that the mass terms can be rewritten: 


1 / 2 


- mlVa) = + rr\V„)(m,<t> - m,Ki) = -4+4- 


( 6 ) 


which dehnes the optimal combinations <h_|_ and <I>_. Since <h_|_ S> for meson models, if 
we rewrite the gradient terms using these variables, we hnd 
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FIG. 4: Neutron skin thickness for three nuclei vs. the isoscalar vector (bottom axis) and scalar 
(top axis) meson masses in a relativistic mean-field model (Ql) from Ref. [^. The circles denote 
the “mean” meson mass (right axis). 
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(top axis) meson masses in a relativistic mean-field model (Ql) from Ref. [^. The circles denote 
the bulk symmetry energy 04 (right axis) in the adjusted parameter set. 
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( 8 ) 


This implies that the “mean” mass fn is much better determined than either of the individual 
meson masses. This conclusion is verihed in Fig. ^ in which m is plotted as circles (the 
scale is on the right axis). The mass m varies by less than 20MeV over the entire range of 
variation of the scalar and vector masses. Thus an (isoscalar) mass scale of about 600 MeV 
underlies the model; we identify it with the scale A used in the power counting. 

One might argue that the usefulness of optimal parameters imply that a non-covariant 
EFT is more appropriate. In Refs. |l^ and this viewpoint is challenged in detail. Here 
we simply note that the large underlying scales characteristic of covariant models are the 
relevant scales in the successful power counting prescription, and that the same scale A is 


consistent with power counting in nonrelativistic Skyrme models . 

If we examine Fig. ^ in more detail, we observe a slight downward trend in the skin 
thickness, particularly for lead; one might think that this is the residual effect of the meson 
range. In fact, we can account for this trend completely by considering another plausible 
candidate for the source of variation in predicted neutron skins: the symmetry energy. We 
introduce the parameter 

N - 7, 

(9) 


a = 


A 


and write the semi-empirical mass formula for the energy per particle in nuclear matter at 
equilibrium as 

E jA = —-|- 040 :^ -f- ttg/A}^'^ . . (10) 

In Fig. the symmetry energy 04 for each reht model is plotted as a circle (the scale is on 
the right axis). The trend in symmetry energy is perfectly correlated with the skin thickness. 

Can most of the variation in neutron skin thickness between models be understood in 
terms of bulk nuclear matter properties, or is the explicit surface symmetry energy also 
important? In Ref. |^, this issue was investigated by considering a macroscopic “droplet- 
type” model of nuclei that divides the energy functional for the binding energy into an 
inhnite matter piece and the lowest-order gradient terms. Specihcally, the binding energy 
R(Z, A) for Z protons and A — Z neutrons was calculated as 

— B{Z,A) = J(frE{p,a)p{r) + Jd^r Eq{\Vp\'^ — PlVpsl"^) + coulomb , (11) 

where p = Pnij") -|- Pp{r) is the baryon density, ps = Pniz) — Pp{r) and a(r) = Pziz)/pij) 
is the local ratio. Surface properties such as the neutron skin come from both density 
dependence in E{p^ a) and explicit contributions proportional to Fq and (3. By htting this 
model to nuclear binding energies, the authors of Ref. concluded that the explicit surface 
symmetry energy (i.e., contributions proportional to (3) was of minor importance. 

This model is quite compatible with the density functional approach. In fact, we have 
precisely the form of Eq. (PH]) with Skyrme and covariant point-coupling energy functionals, 
and for meson models one can use the meson equations iteratively to get estimates for 
the parameters Eq and (3. We can estimate (and bound) the contribution from the surface 
symmetry energy by adopting natural values of the associated parameters and then observing 
the variation in — rp as (3 runs over the natural range. The NDA estimate, which applies 
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FIG. 6 ; Neutron skin thickness in ^osp]^ .^g_ gynimetry energy 04 for a single covariant meson 
mean-field model. The symmetry energy is varied by changing the rho meson coupling constant, 
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FIG. 8: Neutron skin thickness vs. nuclear matter binding energy for a wide variety of mean-field 
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FIG. 11: Neutron skin thickness in ^*^®Pb vs. linear density dependence of the symmetry energy pQ 
for a wide variety of mean-field models, as described in the text and in the caption to Fig. 
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directly to either Skyrme or covariant point-coupling mean-field models, is 


Fo(|Vpr-/3|Vp3r)~/'A' 


LA^ my 


r 


(12) 


Numerical estimates for A = 600 MeV are Fq ~ (hc)^//^A^ ~ 75MeV-fm^ and (3 ~ 0(1/4), 
with the signs undetermined. (The factor of 1/4 originates with the association ^ 

|p 3 .) These values are consistent with all of the Skyrme and point-coupling model fits, 
which is another verification of naturalness. 

Varying (3 between —1 and -|-1 leads to a change in the neutron skin thickness of less 
than 0.03 fm for which is small compared to the typical spread between covariant 

and Skyrme models. Thus the explicit surface dependence of the symmetry energy is not a 
significant factor.^ This means that we can focus our attention on nuclear matter properties. 

In general, the energy per particle of asymmetric matter can be expanded about the 
equilibrium density po in n Taylor series in p and a (we follow the notation of Ref. |p]|: 


E{p, a) = E{p, 0) S 2 {p)a^ + Si{p)a^ H- a = 


K 

0 ) = + T^i.P ~ Pof + ■ ■ ■ 

5*2 (p) = 04+^(p—po) +o 2 (p ~ pyy +■ ■ ■ 
Po i-oPo 


N - Z 


A 


(13) 


(14) 


(15) 


which defines the linear density dependence of the symmetry energy, po, and the correction 
to the incompressibility, AKq. These parameters and the isoscalar equilibrium parameters 
tty, Po, and Kq are our candidates for correlations with the neutron skin. The contribution of 
S 4 (p) for 208 p]^ jg small because of the factor and so 5*4 is not constrained in mean- 
field models, making extrapolations to a ~ 1 (neutron matter) quite uncertain (although 
S '2 dominates near po in any case). We note that a study using realistic interactions found 
S' 4 (p) to be unimportant even at higher densities |]60[ . 

One might try to study the correlations quantitatively by varying appropriate parame¬ 
ters one by one within a single model, as illustrated in Fig. |^. The symmetry energy at 
equilibrium in a single relativistic meson mean-field model (set Q1 from Ref. |M) is varied 
by changing the value of the rho coupling constant Pp [see Eq. (P!B|)]. It is clear from the 
figure that the neutron skin is linearly correlated with the symmetry energy. But if instead 
we plot in Fig. the skin thickness versus 04 for a wide range of models, each with a good 
fit to nuclear properties (binding energies, charge radii, spin-orbit splittings), we see strong 
correlation again but with a significantly different slope. The problem is that one cannot 
vary one parameter independently and expect that the fit remains good. The goodness- 
of-fit value (labeled x^) each set is shown in Fig. with a value of 100 indicating an 
acceptable fit. This shows that the quality deteriorates rapidly unless the other parameters 
are allowed to change. A more sophisticated approach to the correlation analysis during the 
fitting process would be very helpful in future investigations. 


§ There is a loophole here: If f3 is actually determined by long-range pion physics, the scale of the gradient 
might be set by ^ A, which would mean (3 could be much larger than unity. 
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In Figs, p, and the correlations of the skin thickness with the symmetric nuclear 
matter {a = 0) binding energy, equilibrium density, and incompressibility are shown. The 
coarse observation is that there are no dramatic correlations between these quantities indi¬ 
vidually and the skin thickness. At a more hne-grained level one can note: 

• The binding energy of Skyrme models is systematically lower than for covariant models, 
although the difference is small and there is signihcant overlap. 

• The equilibrium density of Skyrme models is systematically larger than for covariant 
models. This has long been observed. There is a very rough correlation with skin 
thickness. 

• The incompressibility seems largely uncorrelated with the skin thickness and there are 
no clear patterns between different types of models. 

Overall, we hnd no signature that the physics determining these quantities directly deter¬ 
mines the skin thickness. Secondary correlations are expected simply from the requirement 
that a good £t is always obtained. 

In Figs. 00 . and |T^, we show the correlations between skin thickness and the individual 
parameters describing the symmetry energy 5*2 (p). In each case the correlation is very strong 
and approximately linear. The spread for po, in particular, is remarkably small and there is 
no signihcant separation of the different models: to good approximation all covariant and 
Skyrme models lie on the same line. A clear corollary to the correlations of 04 , po, and AKq 
to the neutron skin thickness is that these parameters are also highly correlated with each 
other. 


V. DISCUSSION 


Qualitatively, the neutron skin is determined from the energy balance sought for N ^ Z 
between the extremes of equal proton and neutron radii but different densities, and equal 
densities but a sizable neutron skin region. Thus, the cost in energy is determined in 
large part by the density dependence of the symmetry energy, S' 2 (p) |^. In addition, the 
Coulomb repulsion favors larger proton radii, which reduces the size of the neutron skin in 
lead approximately 0.1 fm relative to the radii with Coulomb turned off, with only a small 
residual dependence on the symmetry energy. (This result was obtained by comparing radii 
with and without the Coulomb interaction for each model.) A simple density-expansion 
analysis based on these ideas made by Oyamatsu et al. identifies as a key parameter 
the shift in saturation point with the asymmetry, dpsat/da^ = —9 pq/Kq. The correlation of 
the neutron skin with dp^ax/da^ is shown in Fig. O and is strong, although there is greater 
spread than with pq alone. 

An explanation of the strong correlation in mean-held models between the parameters 
of 5*2 (p) has been ohered by Horowitz |^. There is not enough resolution in the binding 
energy systematics of hnite nuclei to hx S' 2 (p) for a range of densities. In fact, as discussed 
earlier, only one quantity is resolved, so one determines the symmetry energy only at an 
average density (p) for nuclei. (For lead we hnd (p) ~ 0.11 fm~^.) Consequently there are 
many combinations of 04 , po, and AKq (and po) that reproduce this value. This conjecture 
is supported by the observation that the spread of 5*2 ((p)) values among the models used in 
the hgures is several times smaller than the spread of S' 2 (po) values, as shown in Fig. |T^. 
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FIG. 14: The standard deviation of S 2 {p) values for all of the mean-field models as a function of p. 


The standard parameterization of energy functionals for covariant meson models leads to 
an enhanced symmetry energy because the kinetic energy piece is large (since M* is reduced 
from M) and the contribution from the p meson mean field is positive definite. Specifically, 
04 is given by 


1 kl 

6(fc| + M*2)V2 


+ 



Po 5 


( 16 ) 


where Pp is the coupling, nip is the p mass parameter, Pp is the Fermi momentum, and M* 
is the effective nucleon mass ||^. In addition, pq scales with 04 (e.g., the p contribution to 


Po is PpPl/Sm'^p). If fits to nuclei determine the symmetry energy at some averaged density 
less than equilibrium density 5*2 ((p)), this implies that both 04 and po will be large (i.e., a 
large density dependence). In fact, typical values for po are around 6 MeV/fm^, compared 
to about 2MeV/fm^ for most standard Skyrme models, while the corresponding values for 
04 are around 36 MeV and 30 MeV. 

Given this association, we can ask: what should po be? Engvik et ah compared calcula¬ 
tions of asymmetric matter using the best-available nonrelativistic potentials and many-body 
calculations . From the slopes of the curves given in Ref. , one can determine po (up 
to a factor of Pq). One hnds po ~ 2MeV/fm^ with only about a 10% spread. [Note: there 
is some ambiguity in the choice of pg to compare with mean-held models; if a larger value 
associated with some Skyrme models was used, then po would be 10-30% higher.] Two 
Dirac-Brueckner studies of asymmetric nuclear matter |^, ^ yield somewhat higher num¬ 
bers, Po ~ 3MeV/fm^. However, both sets of values are signihcantly less than those from 
the standard covariant mean-held meson models (i.e., the solid black squares in Fig. 0 . 
Generalized meson models can be forced to have lower po (by adding a low neutron radius to 
the ht observables) but sets with pg lower than 5 MeV/fm^ are only achieved if the additional 
isovector parameter (pp from Ref. JT^) becomes highly unnatural. 

Generalized point-coupling models which feature higher-order terms not included 
in conventional covariant models, were found with lower values of po and good hts (see the 
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shaded triangles in Fig. 0). However, this lower value arises from cancellations between 
two orders in the density expansion for the isovector contribution, which again violates the 
principle of naturalness. In particular, the contribution to po analogous to the p meson 
contribution is proportional to C 2 + fhpPO) where po ~ 1/5. By choosing rjp ^ —10/3, one 
can cancel the leading contribution for natural ('2 ~ 1- In fact, one of the sets has Pp = 
—3.245. We do not have a physical reason for such a fine tuning between different orders 
and conjecture that it may reflect instead limitations in the isovector contribution to the 
energy functional. 

The conclusion is therefore that we need to revisit isovector physics in covariant models. 
There are many possible deficiencies to explore: 


The functional is incomplete. For example, isovector scalar contributions (the 
meson) may be important for isovector properties. This possibility has been recently 
considered in Ref. . 


• The functional is truncated too severely. The higher-order terms discussed above are 
an example of what could be added. 


• Long-range physics must be explicitly included, as expected in an EFT treatment. 
One natural class of contributions is from the pion. As noted above, explicit pions 
have not been necessary because the isoscalar part is shorter-ranged (“a” physics) 
and the isovector part is poorly resolved. However, power counting implies that pion 
contributions should be included. In addition, many-body long-range correlations 
(e.g., RPA correlations) may be needed. 

• The one-loop form of the functional is not sufficiently general. For example, one will 
obtain nonanalytic terms in the density from the long-range contributions. 


We seek a systematic approach to incorporating this physics into our descriptions of medium 
and heavy nuclei. Furthermore, a connection to underlying forces is desirable to help con¬ 
strain the functional. 

A key question is whether mean-field models can be made more systematic without dras¬ 
tically increasing the computational burden. A great appeal of mean-held models is the ease 
with which they can be used to calculate medium to heavy nuclei. The basic calculational 
procedure is the same for a covariant Hartree or nonrelativistic Skyrme interaction. The 
models can all be formulated in terms of an energy functional, whose minimization with re¬ 
spect to single-particle orbitals yields a single-particle “potential” functional. This potential 
is used in turn to calculate the orbitals. 

The generic procedure is quite simple. To solve for the ground-state bulk properties of a 
nucleus with A = Z + N nucleons. 


1. guess a set of initial density/current prohles {p(x)} (this notation denotes the baryon 
and isospin densities, plus any other spin or scalar densities); most can be set to zero 
or initialized with a Fermi shape with appropriate parameters; 

2. evaluate a functional of the {p(x)}, yielding a local single-particle potential IA(x); 

3. solve the Dirac or Schrodinger equation for the lowest A eigenvalues and eigenfunctions 

[-ia. ■ V + f3M + f3Vs{-x)] = eai>a{^) , (17) 
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r V2 
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+ (x) ^pa{x) 


e«V'Q(x) , 


4. compute new densities, e.g., 


(18) 


A 

q ;=1 


(19) 


and other observables (which are functionals of {ea, V'a}); 

5. repeat steps 2.-4. until the changes from iteration to iteration are acceptably small, 
i.e., until the solution is self-consistent. 


This procedure is straightforward to implement numerically. Dealing with deformed 
nuclei and open shells introduces some complications (e.g., pairing), but these do not qual¬ 
itatively complicate the calculations. There are two key reasons for the simplicity: 


• the energy functional is universal, in the sense that the same functional is used for all 
nuclei, with the same set of parameters; 

• the single-particle functional, which is derived from the functional derivative of the 
energy functional with respect to p(x) is local (at least when pairing is treated most 
simply) and is easily evaluated in terms of the eigenvalues and eigenfunctions. 


Once the functional is established, each subsequent hnite nucleus calculation is almost trivial. 
Thus it is certainly desirable to maintain this mean-held calculational procedure. 

In Ref. P] , it is argued that the mean-held approximation for meson-ranged interactions 
cannot be valid at ordinary densities, because the inverse Compton wavelengths /i of the 
mesons times the mean interparticle spacing is much less than unity. The authors go on 
to state that “The RMF approximation can be based on ehective values of the coupling 
constants that take into account the correlation ehects. However, these coupling constants 
then have a density dependence, and a microscopic theory is needed to calculate them.” The 
second objection is misleading, since density dependent couplings are not necessary; all of 
the density dependence can be incorporated in the functional, if made sufficiently general. 
The correspondence between models with explicitly density dependent couplings and more 
general functionals can be made through held redehnitions [^. Matching to a microscopic 
theory may be a desirable way to determine the parameters, but it is not necessary. Instead, 
one can match to hnite density observables, which is especially needed for calculating heavy 
nuclei. 

To address the hrst objection, we propose a merger of Density Functional Theory (DFT) 
and Ehective Field Theory (EFT). Kohn-Sham DFT is a framework that looks just like 
the simple relativistic Hartree calculations that are so well adapted to calculating hnite 
systems, but one which can include all short-range and long-range correlations |^. That 
is, the solution framework described above can accommodate the most general description 
of the ground state within the context of density functional theory, and the Kohn-Sham 
scheme in particular. The implication is that conventional mean-held approaches provide 
reliable descriptions of bulk properties not because the mean-held approximation is good but 
because they accurately approximate exact density functionals. This claim is supported by 
the work of Hu, who showed that Dirac-Brueckner-Hartree-Fock at the two hole-line level. 
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when matched to a generalized EFT-based functional, is reproduced with natural coefficients 

H- 

Therefore, we interpret the “mean-held” functionals as just particular generalized gradi¬ 
ent approximations to complete Kohn-Sham density functionals ||I3|, [^. (Skyrme may be 


considered more general in this regard.) Rather than starting from the underlying free-space 
interaction, one expands the functional as a series expansion in density and momentum, and 
then hts directly to experimental data. This approach can miss (long-range) non-analytic 
or nonlocal contributions, such as explicit exchange (“Fock”) terms (e.g., from pions). The 
role of the EFT framework is to add these systematically while maintaining symmetries and 
conservation laws. 

From a Green’s function perspective, the mean-held models follow from the Hartree (rel¬ 
ativistic) or Hartree-Fock (nonrelativistic) approximation to the nucleon self-energy, with 
an ehective interaction (e.g., an approximation to the G-matrix). In practice, the non¬ 
relativistic Skyrme interaction uses zero-range (“point couplings”), which means that the 
energy functional has the same form as for a Hartree calculation (one can make a Fierz 
rearrangement of the interaction to derive the correspondence). The Hartree self-energy in 
coordinate-space is static and local (i.e., a function of x only). 

In contrast, the general self-energy is non-local in space and is frequency dependent. This 
would seem to imply that the correlation ehects that go beyond the Hartree approximation 
could never be systematically included in the calculational framework described above; one 
could only incorporate averaged ehects using ehective interactions. Since the simplicity of 
the mean-held approach depends heavily on the local, static nature of the single-particle 
potential, it would appear that systematic improvement is not possible. The answer is that 
we don’t need to calculate the general self-energy. Instead we use the Kohn-Sham DFT 
potential 14 (x), which is static and local, even though it can include all correlations. 

There is a danger that evaluating the single-particle functional to hnd 14 (x) will become 
too great a computational burden when going beyond the current mean-held models. So 
developing reliable gradient expansions will be essential. The surprising experience from 
Goulomb systems that LDA plus generalized gradient expansions work extremely well needs 
to be understood in this context. In particular, it does not rule out long-range physics being 
incorporated in the same form as existing functionals, although contributions nonanalytic in 
the densities should be expected. For example, the exchange contribution from a long-range 
(i.e., ehectively massless) particle will have a contribution to the energy density. Thus, 
an EFT treatment of the pion expanded around the chiral limit should contribute a term 
like this. Work is in progress to cast Kohn-Sham DFT into an ehective action formalism for 
composite operators . 


We could imagine constructing covariant models more like Skyrme models, by including 
terms like pf and p^, with fractional S and 7 determined phenomenologically. But the 
construction from a Lagrangian has been an important guide to maintaining covariance and 
building a conserving approximation. In addition, the connection to free space scattering 
is much closer than in the nonrelativistic case since self-energies are dominated by mean- 
held components (Hartree dominance), so corrections are smaller. Therefore, we advocate 
deriving a DFT expansion from an EFT Lagrangian, or at least to establish constraints on 
the analytic structure. 

Finally, we return to the issue of point nuclear densities. Single-nucleon structure will 
automatically be included in energy functionals ht to data. Since point proton and neutron 
radii are not observables (no probe couples to them) they are at best auxiliary quantities. 
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This is not a major issue for properties such as the charge density alone, since the nuclear 
structure is mostly isoscalar while the single-nucleon form factors are mostly isovector and 
so factorization is a good starting point. But considering — rp heightens sensitivity to the 
isovector form factor. In the future, one should actually couple to electroweak currents to 
avoid the model dependence inherent in defining the point densities. One of the strengths of 
an EFT is knowing how to couple to external currents, which may involve identifying new 
constants (low-energy constants). 


VI. SUMMARY 


We have examined the predictions by mean-field models of neutron skin thicknesses, 
which are highly correlated with the neutron radius. We find that the skin thickness is 
closely related to the density dependence of the symmetry energy, independent of the type 
of mean field model. In particular, the linear density coefficient po exhibits a remarkably 
good linear correspondence. (However, the strong correlations between parameters preclude 
isolating the effects of a single parameter.) This implies that it is a general feature of the 
energy functional. Our conclusion is in accord with other analyses in the recent literature 


m 


An accurate extraction of the radius from experiment, for even a single nucleus, would 
provide a new and valuable constraint on the energy functional. Current analyses of proton 
scattering are consistent with a skin thickness in lead of order 0.1-0.2 fm which would 


imply 1 MeV/fm < po < 4 MeV/fm , and recent analyses of antiprotonic atom data also lie 
in that range . Calculations of uniform matter with realistic potentials show a narrower 
range of 2 MeV/fm^ < po < 3 MeV/fm^, with Dirac-Brueckner calculations at the high end 
and nonrelativistic calculations at the low end. If these conclusions are validated, it means 
that the uncertainty in neutron radii may be much smaller than implied by comparisons of 
the current “best fit” calculations. 

However, the commonly used mean-field functionals are not equally flexible in accommo¬ 
dating different values of Pq. In particular, covariant meson models have significantly larger 
values of po than found in most Skyrme parametrizations, such as the recent set of Brown 
Ij^. Calculations of uniform matter with realistic potentials imply that po is close to that 
of nonrelativistic Skyrme models while it is significantly overestimated in the most widely 
used relativistic mean-field models. It is not surprising that this does not preclude a good 
fit to the best measured properties of hnite nuclei, given the limited constraint from these 
properties on purely isovector parameters. It is not clear whether the tendency for Skyrme 
models to have lower po is a consequence of the form of truncation or is slightly favored by 
the fits to other observables. In Ref. additional information, which corresponds to fixing 
Po, was used to constrain one of the isovector parameters. 

Such comparisons between nonrelativistic and covariant models might be used to argue 
that one approach is more “correct” than another. However, the EFT perspective changes 
the context of the discussion from correctness to efficiency and the effects of truncation. 
The most commonly used covariant parametrization appears to be deficient because it does 
not allow lower po with natural parameters, but the conclusion should be that the func¬ 
tional must be improved. Power counting suggests several avenues: a more complete set 
of operators at the lower orders (e.g., add an isovector scalar meson or density), higher 
order contributions (density corrections to rho meson exchange), pion and other long-range 
contributions. Additional isovector observables, such as binding energies of more asym- 
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metric nuclei and observables from collective excitations will be needed to constrain new 
parameters. The question is: how can we incorporate this physics to improve the covariant 
functionals (and nonrelativistic functionals) systematically? We propose developing density 
functional theory in an EFT framework, using effective action techniques |p5|| . 
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